Benedicks effect in a relativistic simple fluid 

A. L. Garcia- Perciante^, A. Sandoval- Villalbazo^, L. S. Garcia-Colin^ 
^Depto. de Matematicas Aplicadas y Sistemas, 
Universidad Autonoma Metropolitana-Cuajimalpa, 
Artificios 40 Mexico D.F 01120, Mexico. 
"^Depto. de Fisica y Matematicas, Universidad Iberoamericana, 
Prolongacion Paseo de la Reforma 880, Mexico D. F. 01219, Mexico, and 

^Depto. de Fisica, 

^ I Universidad Autonoma Metropolitana-Iztapalapa, 

Av. Purisima y Michoacan S/N, 



O 

^ : Mexico D. F. 09340, 



Mexico. Also at El Colegio Nacional, 
Luis Gonzalez Obregon 23, 



O ■ Centro Historico, 

^ . Mexico D. F. 06020, Mexico. 

> . Abstract 
(N 

^ : According to standard tlrermophys:cal theories, cross effects are mostly present in mnlt:compo- 

' nent systems. In this paper we show that for relativistic fluids an electric field generates a heat flux 

(N 

' even in the single component case. In the non-relativistic limit the effect vanishes and Fourier's law 

m 

is recovered. This result is novel and may have applications in the transport properties of very hot 

; plasmas. 
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I. INTRODUCTION 



It is well known that when one considers an ionized plasma as a binary mixture in the 
presence of an electromagnetic field, direct and cross effects arise. Of special interest, are 
the thermoelectric effects which have been discussed in Ref. |]|. The Benedicks effect, also 
called electrothermic effect, corresponds to a heat flux caused by an electric fleld. It's history 
can be traced back to the experimental works developed by Carl Axel Fredrik Benedicks in 
the second decade of the twentieth Century in homogeneous systems [2]. These experiments 
showed the existence of the reciprocal cross effect of the Thomson (thermoelectric) effect 
(1856) which consists of the generation of a flow of charge due to a temperature gradient. 
The Benedicks effect should not be confused with the Peltier effect, which is strictly related 
to inhomogeneous systems. 

Here, we present what we consider an interesting result which may lead to applications. 
If one takes a single component charged relativistic gas subject to the presence of an electric 
field, in addition to the electric current a heat fiux appears driven by such field. This effect is 
absent in the non-relativistic domain and suggests the possibility of having additional cross- 
like effects in relativistic fiuids. This effect may also impact the dynamics of hot plasmas 
possibly affecting stability and other properties of such systems. 

Several authors have addressed the problem of heat conduction in charged fluids. In Ref. 
jsl a calculation of the thermoelectric effect, together with the electric current associated 
with a temperature gradient is performed using Anderson and Witting's collisional model 
for a binary mixture. Also, in the non-relativistic scenario a reciprocal effect giving the 
electrical current due to a temperature gradient is usually obtained assuming stationary 
state in Boltzmann's equation [4], a procedure that will be carefully discussed in the last 
section of the present work. 

In this work we establish the relativistic heat flux for a simple charged fluid in the presence 
of an electric fleld starting from a relativistic Boltzmann equation using Marie's relaxation 
time approach to describe the collisional effects. To accomplish this task, a correction to 
the relativistic local equilibrium distribution function is obtained in Sect. II in terms of the 
gradients in the system to flrst order in the Knudsen parameter. Section III is devoted to 
the explicit calculation of the heat flux and the corresponding transport coefficients. The 
interpretation and discussion of these results are given in Section IV together with final 
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remarks. 



II. THE RELATIVISTIC BOLTZMANN EQUATION 



We start the calculation by setting up the kinetic equation for the system in consideration 
which is a dilute gas of non-degenerate charged particles at a temperature high enough so 
that z = is at least close to unity and relativistic effects in the individual dynamics are 
relevant. The system as a whole has a hydrodynamic four velocity which we denote by 
and its constituents have charge q and mass m. The relativistic Boltzmann equation in the 
relaxation time approximation reads 
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The structure of the BGK-like collision kernel was first suggested by Marie and is consid- 
ered here as an adequate approximation since it will lead to the general structure of the 
distribution function as well as the dissipative fluxes without introducing all the compli- 
cated expressions and calculations that are involved when considering the complete collision 
kernel. The complete calculation would perhaps yield relevant corrections to the particular 
values of the transport coefficients but both the structure of the fluxes as well as the general 
behavior of the coefficients will remain unaltered. Here / is the distribution function, f^^^ the 
local equilibrium solution, and r a characteristic relaxation time. The molecular dynamic 
variables are the individual particle four-velocity and the corresponding acceleration v^^, 
both referred to the laboratory frame. 

Since the fluid is charged, in presence of an electrostatic field the acceleration in the 
second term on the left side of Eq. ([T]) may be written as 
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is the usual electromagnetic field tensor and thus 

= 1,2,3 

mv^ = ^ ^ ' (4) 

/i = 4 

Here a coma indicates a partial derivative whereas a semicolon will denote a covariant one. 
For the special relativistic fluid, we consider a flat Minkowski spacetime with a + + H — 
signature such that the covariant derivative will coincide with the ordinary one in cartesian 
coordinates. The dot represents a proper time derivative. 
Substitution of Eq. (j2]) in Eq. ([1]) leads to 

+ -VaF-'^^ = (5) 

which will be solved, following the standard Chapman-Enskog procedure, by introducing 
the assumption that the distribution function can be written as 

/ = /(°^ + /^^) (6) 

where Z*-^-* is of flrst order in the gradients. By substituting Eq. (jH]) in Eq. ^ and keeping 
only terms up to flrst order in the gradients, one can write the deviation from the equilibrium 
solution as follows 

f,..r[n-'^{^^.F''Y-g,] (7) 

The presence of the local equilibrium distribution function on the right side of Eq. ([7j) 
follows from the fact that its derivatives are already flrst order in the gradients since the 
local equilibrium assumption asserts that its space and time dependence is given only though 
the gradients and time derivatives of the state variables n, T and W . It must also be taken 
into account that the time derivatives of such variables are in turn written in terms of 
gradients by introducing Euler's relativistic equations, a step necessary in order to assure 
existence of the flrst order solution jo, Q|. 

In order to establish the explicit dependence of /"^^^ with the gradients, we start by using 
the following identity 

t^" = t;^/^^ + 7wW° (8) 

where h'^ projects in the hyperplane orthogonal to the hydrodynamic velocity, that is h'^lA^ = 
0. Also, for the last equality use has been made of the identity f'^W/? = — 7(fc)C^ which is 



valid since, being the contraction a scalar, it can be calculated in a comoving frame where 



0,c 



and = K"", with 



ir" = 7w[A;^c] £ = 1,2,3 (9) 

denoting the chaotic or peculiar velocity and 7(fc) = (1 - k'^/c^) being the usual 

Lotentz factor with k = \fWkii. By noticing that f is also an invariant quantity, the 
calculation is carried out in the comoving frame where the projector has the simple form 

/l o\ 
10 
10 
\0 0/ 



(10) 



so that one can write in general 
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where i, i = 1,2,3 and the derivatives are calculated in the comoving frame, denoted by the 
subscript CF, where local equilibrium is satisfied. Since locally the equilibrium distribution 
function is a Juttner function given by 5|, l8| 



f(0) = !1 exp^^^ 



one has 
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where Q (i) = ^ and /C„ (i) is the modified n-th order Bessel function of the second 
kind depending on the relativistic parameter z = . Thus, 
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Following Hubert's procedure the time derivatives are to be written in terms of the 
gradients of the state variables via Euler's (lower order) equations: 



(17) 
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where C„ is the specific heat for constant n and for and ideal gas p = ^ + ^ = nmQ (i). 
The set given by Eqs. ffT7l) - ffT^ corresponds to the local equilibrium, zero order in the 
gradients, relativistic hydrodynamic equations which can be obtained in the kinetic theory 
framework by multiplying Boltzmann's equation by collisional invariants and integrating 
in velocity space. Such a procedure can be found for the relativistic case in Ref. js] for 
example. In order to obtain such relations in the comoving frame, as required by Eq. (fT6l) . 



we consider W 
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and Eq. (fTO|) which leads to 
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Introducing Eqs. fl2U|) to f l25]) in Eq. (ITB]) one obtains 
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where a term proportional to the velocity gradient is ignored from now on since it will not 
couple with the heat flux. This fact is sustained by Curie's principle by means of which only 
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fluxes and forces of the same tensorial rank are coupled to each other in the constitutive 
equations. 

On the other hand, the second term in Eq. ([7j) is calculated using the fact that the fleld 
tensor only has electric potential components (4 — £ with ^ running up to 3) such that 

^..F-) ^ = = ^..F-^ (25) 

m / zc^ m zc^ m zc^ 

The quantity is scalar and therefore can be calculated in the comoving frame, so that 



-m / zc'^ m zc 

and thus, since = i— = _ 



^ = 7(.)^fc.F^^— (26) 




, / 1 7(fc) l{k) Q (i) \ T, [ 

The function in Eq. (127|) is the correction to the equilibrium distribution function to flrst 
order in the gradients in the presence of an electric fleld. This part of the distribution yields 
the dissipative terms in the hydrodynamic equations to the Navier-Stokes level. In the next 
section, f^^^ will be introduced in the expression for the heat flux in the relativistic case in 
order to establish the general structure of the corresponding constitutive equation as well 
as the particular value of the transport coefficients for the relaxation approximation. 



III. CONSTITUTIVE EQUATION FOR THE HEAT FLUX 

An analysis of both the concept and the formal expression for the heat flux in relativistic 
cinetic theory in relation to the flux of chaotic energy, as initially conceived by R. Clausius 
qI, can be found in Ref. [1]. In this context it can be shown that the heat flux in the 
comoving frame is given by 

JlQ] = nicX I k^f^hl)<rK (28) 

and can be expressed in an arbitrary frame using a Lorentz transformation. From fl28|) and 
considering the structure of Eq. (127|) it can be seen that three independent driving forces 
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will be present in the flux namely, the temperature, particle density and electric potential 
gradients. Of particular interest in this work is the term that corresponds to the electric 
fleld which in this relaxation time approximation is given by 

Jkm = -r'^ j - 1) o4,<;-/f. (29) 

Notice that this flux vanishes in the non-relativistic limit where both 7^ and Q (i) are unity. 
Thus, the heat flux in the single component fluid due to the thermoelectric effect is only 
present for non-negligible values of the parameter z and can be written as 

Jf^,^] = LteKV (30) 
where the corresponding transport coefficient is given by 

(31) 

\-z) j 

For small z, in the non-relativistic limit one has 

Lte = riTc^q^ {z" - z"^ + ...) (32) 

or equivalently 

I (^^^ 
The general form for the heat flux including all the gradients is given by 

Jfg] = -h^"" (ltt^ + Ltu^ + LTE(t>,^ (34) 
where the transport coefficients corresponding to the first two terms were previously obtained 




12|. In the non-relativistic limit, the ordering of the effects with respect to c can be extracted 



from the following expression 
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T ' mc^ n mc^ ' 

from where it can be clearly seen that the only effect that survives in the non-relativistic 
limit is the one corresponding to Fourier's law and the other two terms are relevant only 
when the thermal energy and/or the electric energy are comparable in magnitude to the rest 
mass of the molecules. It is worthwhile to remind the reader at this point that the results 
obtained in this section are valid for single component gases, no binary mixture needs to be 
considered. 
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IV. SUMMARY AND FINAL REMARKS 



In this paper, we have studied the heat flux in a single component relativistic charged 
gas in the presence of an electric fleld. The calculation was performed using Marie's kernel 
for the Boltzmann equation. The result is a heat flux coupled to three forces, namely the 
gradients of the temperature, density and electric potential. Only the temperature gradient 
contribution does not vanish in the non-relativistic limit while the other two effects are 
strictly relativistic. 

Some authors have argued in favor of a cross-like relationship between temperature gra- 
dients and electric fleld in the non-relativistic single component fluid [4|. In order for such 
an effect to be established, one would have to set the partial time derivative in Boltzmann's 
equation to zero. This is usually justifled by assuming a steady state in the gas, however 
the Chapman-Enskog procedure arises from an ordering scheme in Boltzmann's equation in 
which the time derivative of the equilibrium distribution function results as a flrst order in 
the gradients term and thus must be considered when calculating the flrst order correction. 
Also, if the Benedicks effect would exist in the non-relativistic single component fluid in 
steady state, one should expect it to be recoverable as a limit of the complete solution. 

By inspection of the structure of the heat flux given by Eq. (134|) . one can draw inter- 
esting conclusions regarding the physics of dissipative simple relativistic charged gases. If 
for the sake of simplicity one ignores the particle density for a homogeneous system, the 
constitutive equation has a structure similar to the one leading to Tolman's law in the case 
of a gravitational fleld [13]. Reasoning in a similar way as in Ref. [13], if one assumes a 
vanishing heat flux, then to lowest order in z: 

^ = (36) 

in such a way that a temperature gradient can be sustained in the absence of a heat flux 
by means of an electric fleld for the simple charged gas. It is interesting to notice that for 
the case of negative and positive charged particles, the required electric potential gradient 
reverses direction for the same temperature gradient. 

Another interesting question arises here regarding the possibility of a relativistic version 
of a Wiedemann-Franz like relation for a single component fluid. This could lead to exper- 
imental tests for the strictly relativistic effects here presented, as well as for the pertinent 
generalizations of this formalism including magnetic flelds. 
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